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We present an analytic derivation of the linear temperature dependence of the inverse static 
susceptibility x~^{T, U) ^ T ~ Tc{U) near the transition from a paramagnetic to a ferromagnetic 
correlated metal within the dynamical mean-field theory (DMFT) for the Hubbard model. The 
equations for the critical temperature and interaction strength of the transition are also determined. 
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I. INTRODUCTION 

The dynamical mean-field theory (DMFT) is a non- 
perturbative and thermodynamically consistent approx- 
imation scheme for quantum mechanical many-body 
problems on a lattice (for reviews see Ref . ^, || ) which 
becomes exact in the limit of large coordination numbers 
[Q. In contrast to the static Hartree-Fock mean-field 
theory the dynamics of the quantum mechanical corre- 
lation problem is fully included in the DMFT. In the last 
few years the DMFT proved to be a powerful tool for 
the investigation of fermionic lattice models with local 
Coulomb interaction such as the Hubbard model and the 
periodic Anderson model [0, |^. It is particularly use- 
ful in the case of intermediate-coupling problems such 
as the Mott-Hubbard metal-insulator-transition or itin- 
erant ferromagnetism where perturbative techniques fail 
(for reviews see Refs. fl], |, §)■ In the DMFT the lattice 
problem is mapped onto an effective single-site problem 
whose self-energy S(ci;) and Green function G{uj) have 
to be calculated self-consistently with the k-integrated 
Dyson equation. The theory is therefore purely local, i.e., 
the self-energy S(ijj) is k-independent and the propagator 
Gk(w) = — ^{^)) niay be represented by the non- 

interacting propagator at shifted frequency, at least 
in the paramagnetic case. Here the mean-field charac- 
ter of the theory becomes particularly evident. The local 
nature of the theory implies that short-range order in 
position space is missing. 

Numerical solutions of the DMFT equations, in par- 
ticular by quantum Monte-Carlo simulations (QMC), re- 
vealed that in the case of continuous phase transitions 
(e.g., from a paramagnetic metal to an antiferromagnetic 
insulator or to a ferromagnetic metal Q) the static 
susceptibility x{T) shows a Curie- Weiss behavior above 
Tc, i.e., oc T - Tc, implying that x oc (T-Tc)"^ 
diverges with a critical exponent 7 = 1. Furthermore, 
the order parameter was found to vanish with an expo- 
nent /3 = 1/2 ^. In view of the mean-field nature 



of the DMFT these numerical findings did not come as 
a surprise. However, considering the dynamics of the 
quantum mechanical problem the situation is not as self- 
evident as it may seem. At a continuous phase transition 
between two phases in high dimensions, a mean-field be- 
havior is naturally expected if the low energy excitation 
spectrum of both phases has a gap. In this case, the 
fermions may be integrated out, leading to an effective 
Ginzburg-Landau-Wilson field theory This applies 
to transitions in the Heisenberg spin model. By contrast, 
if the transition occurs between two metallic phases, e.g., 
from a paramagnetic to a ferromagnetic metal, or be- 
tween a metallic and an insulating phase, the result is far 
from trivial since the low-lying excitations in the metallic 
phases may couple to the order parameter and thereby 
lead to divergences in the effective field theory at T = 
or even at very low temperatures ||ll|] . At present, the 
consequences of this feedback are still not entirely under- 
stood. It is therefore worthwhile to further investigate 
the transition between two metallic phases, also within 
the DMFT, using analytical methods. 

Curie-Weiss behavior of the magnetic susceptibility 
is traditionally associated with localized magnetic mo- 
ments, and, indeed, is readily obtained for Heisenberg- 
type spin models in a mean-field approximations. Nev- 
ertheless, it is also a characteristic of interacting itiner- 
ant electrons as described, for example, by the Hubbard 
model. In particular, a Curie- Weiss behavior may be 
obtained within the Hartree-Fock approximation which 
yields X^U^) ~ - above § |l|. However, 
since this result is derived for interaction strengths where 
the Hartree-Fock approximation is not controlled by a 
perturbation theory, its qualitative and quantitative va- 
lidity is questionable The same criticism applies to 
the Stoner criterion for the onset of ferromagnetism . 

To calculate the static magnetic susceptibility in the 
vicinity of a continuous phase transition from the para- 
magnetic to the ferromagnetic metallic state where elec- 
tronic correlations are explicitly included we will use the 
one-band Hubbard model 

^ ^Yl ^ycL^j^ + UY1 '^iT'^U' (1) 

ij,(7 i 
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on an arbitrary lattice and employ the DMFT. In par- 
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ticular, we will show that the critical exponent is indeed 
7 = 1. 



II. DERIVATION OF THE SUSCEPTIBILITY 

We wish to calculate the magnetization density 

TO = (1/2/3) ^aG,„, (2) 

n 

where the local Green function Gan in DMFT is given by 
the bare density of states N^{() and the local self-energy 



Go 



de- 



iujn + IJ, + ah — J^o 



(3) 



Here the subscript n refers to the Matsubara frequency 
itOn = i{2n + l)7r//9 for the temperature T, with /3 = 
l/ksT, and h is the external magnetic field in energy 
units. Within the DMFT the local Green function Gem 
is determined self-consistently through (|^) and 
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by the k-integrated Dyson equation 
The single-site action A has the form 

n.cr 

U / dTcUT)cAr)4ir)c^{T), 
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(4) 
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where we used a mixed time/frequency convention for 
Grassman variables Co- , c*. 

We first separate the self-energy So-n into its static (i.e., 
Hartree-Fock) part and its explicitly dynamical contribu- 
tion Eo-n as 



San = t^-TT — all TO + 



(7) 



Here uq is the density of particles. The Hartree-Fock 
approximation corresponds to neglecting Eo-n- Since we 
are interested in the behavior of the susceptibility close 
to a continuous transition, i.e., in the limits T — + Tc, 
/i — + 0, TO — > 0, where is the, yet unknown. Curie 
temperature, we write 
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Gn + SG 

an : 
Gn + 5Qan, 
E„ -|- SY^om 



(8) 



where SGan, SQan, and S'Ecn are spin-dependent cor- 
rections to the respective paramagnetic (i.e., spin- 
independent) parts. In the next step we expand (^ and 
(^) up to first order in these corrections. Using Eqs. (||) 
and (^) we find 



SGan = - [(Th + aUm - 5Y.„i]Fn, 
where we introduced the function 

iV"(e) 



Fn 



de- 



(9) 



(10) 



As in the Weiss molecular field theory, the particles may 
be interpreted as moving in an effective magnetic field 



heff — h + Um — crSHa 



(11) 



However, in the case of interacting electrons considered 
here this effective field is, in general, found to be dy- 
namic, i.e., it fluctuates in time due to the local cor- 
relations which lead to an exchange of energy between 
the particles. Neglecting the dynamical term (5Ecrn we 
recover the usual static mean-field expression for the ef- 
fective magnetic field. Similarly, expanding Eq. (^ and 
using Eq. (0) we obtain 
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e-^ = G-^ + E„ + C/^, 



and 



SGan SG 

an 



Gn 



aUm — (5En 



(12) 



(13) 



Our goal is to find S'Ean and then, with the help of En . 
(^), to calculate the magnetization to. Since Eq. (I14 ) 
still contains the two unknown functions STi^n and SQa-m 
an additional condition is required to close the set of 
equations. This condition is provided by the fact that the 
self-energy is a functional of Q^n in perturbation theory, 
i.e., E(jn = 5jcr[^(Tn], to infinite order. Hence, we can 
formally expand as 

Scr [Ga7i] = Scr [Gn + SQan] « E [Qn] + J^"^ SQ^n' ■ 

, oyn' 

n' 

(14) 

Using Eq.(|l^ relating Q and G we find the functional 
derivative 



(5S _ ^E 5G I 5G 

'sg ' SG "sg ^ F'sg' 



(15) 



where F = F„ is given by Eq. (^. Employing Eq. (Q) 
in the paramagnetic phase and using Eq. (p^, we obtain 
the spin-dependent correction to the self-energy Ecm as 
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where 



and 



-a [/ii/„ + C/m(i/„ + Ai/„)] F„ 



/ J rin' / J TP 



E-"- n'n' 



hHyi + UmiHn + AHn). 



Gl, 



(16) 



where M^^, is the inverse of the matrix Af„„' defined as 
M„„, = 5„„, + _ 1^ , (17) 

and 



is the two-particle density-density correlation function 
calculated in the paramagnetic phase. We note that in 
the Hartree-Fock approximation the two-particle correla- 
tions are neglected, i.e. F^^^ = 0, and therefore So-n = 0. 

The spin-dependent correction to the local Green func- 
tion Gan can now be expressed as 



(19) 



(20) 



(21) 



We see that the effective magnetic field ([ll| ) acting on an 
electron is given by 



(22) 



It is interesting to observe that the dynamics of the two 
terms is different^ i.e., the correlation problem leads to 
an asymmetry between the external (h) and the induced 
(Urn) effective magnetic fields. The origin of this asym- 
metry lies in the self-consistency equation (^, where 
h enters through SG^n , while U m enters both through 
SGan and S(j„. In the Hartree-Fock approximation, the 
frequency dependent factors iJ„ and A_ff„ reduce to unity 
and zero, respectively, such that ft-e// — h + Um becomes 
a static effective magnetic field. 

We are now able to calculate the magnetic suscepti- 
bility x^'^^^'^ {T,U). Noting that the magnetization m, 
Eq.(2), has nonvanishing contributions only from 5Gnm 
and using Eq.([l9|) one finds 

^-Y,HnFnh-'^Y.^H,, + AH^)FnUm. (23) 
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The linear magnetic susceptibility is then obtained from 

^^^DMFT^rj.^ [/)/iaS 



x^^^^^(r,t/) 



Xo(T, U) 



(24) 



where 



and 



Xo 



l-C/[xo(T, t/)-f Axo(r, U)] 
{T,U) = --Y.HnFn, (25) 



Axo(r, t/) = -i^Ai7„i^„ 



(26) 



Eq. ( p4| ) with Eqs. ( |25| ) and ( p6| ) is one of the main re- 
sults of our paper. We note that the expression for the 
static susceptibility x"^'^^"^ {T,U) in Eq. (|2|) is decep- 
tively simple. Indeed, it has the form of the correspond- 
ing RPA expression, with the Pauli susceptibility Xo{T) 
of the non-interacting system replaced by the susceptibil- 
ities xo{T,U) and Axo{T,U) of the interacting system. 
If both E„ and F„„/ were neglected we would recover the 
well-known Hartree-Fock result. The result for the static 
susceptibility can be expressed in the RPA-like form 



X 



,DMFT 



{T,U)^ 



X{T, U) 
1 - Ux{T, U) 



with 



x(r, u) = 



Xo(T, u) 



i + c/Axo(r, [/)■ 



(27) 



(28) 



This equation suggests that, as in RPA, the suscepti- 
bility can be written as an infinite series of bubble dia- 
grams with (3-leg) vertex corrections. So far it was not 
possible to deconvolute the corresponding Bethe-Salpeter 
equation for the vertex corrections, since the scattering 
function (4- leg vertex), although is k-independent in the 
infinite dimensional theory, is still a complicated func- 
tion of frequency. It should be noted, however, that our 
algebraic derivation of the static susceptibility is non- 
perturbative anyway since it is not based on direct dia- 
grammatic resummations. 



III. DISCUSSION 

The transition point between the paramagnetic and 
ferromagnetic phase is determined by the divergence of 
the static susceptibility x^^^^^(T, U) in Eq. (|2|) 



l-C/[xo(re,C/) + Axo(Tc,C/)] =0. 



(29) 
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Using the spectral representation for the summation over 
the Matsubara frequencies one can write (E9[) in a closed 



J 



form as 



1 - 



de 



1 



Im [H{e + iO+) + Ai7(e + iO+)] \ = 0, 
I 



(30) 



where H{e + iO~^) and AH{e + zO+) are obtained by 
analytic continuation: 77„ = H{iijjn) H{e + iO^) 
and AHn = Ai7(iw„) -> AH{e + iO+). It should be 
noted that H{e) and AH{e) are still functions of tem- 
perature because of the internal summations over Mat- 
subara frequencies in Eqs. (20) and (|2^). For given 
U Eq. (30) determines the Curie temperature Tc{U). 
Similarly, one may fix the temperature to determine the 
critical interaction strength Uc{T). In this respect Eqs. 
(0), ( ^o| ) are generalizations of the Stoner criterion |l|] 
Ijjtoner ^ I/N^Uf) obtained in Hartee-Fock theory In 
general Eq. (^ implies that, due to the inclusion of 
genuine correlations, the transition point to the ferro- 
magnetic phase is not merely determined by the density 
of states at the Fermi level but, rather, by the density of 
states at all energies. Due to the increase in the kinetic 
energy the value of Uc is reduced by an asymmetric DOS, 
especially if the DOS has a singularity at the lower band 
edge. This had already been found in the numerical so- 
lution of the DMFT equations for a model DOS using 
Monte-Carlo simulations Q and in the approximate 
treatment of DMFT within the modified perturbation 
theory Eq. ( pO| ) confirms these numerical findings 
analytically and explains the origin as a correlation effect. 

Expanding the static susceptibility ( p4| ) around Tc{U) 
(for T > Tc > 0) we find a Curie- Weiss law 



,DMFT 



{T,U)^ 
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{Tc{U),U) 



with 
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iTciU),U) 



T - T,{U) 



Xo{Tc{U),U) 



(31) 



jj (^ d\xoiT,U)+Axo{T,U)] j 



(32) 



which is hence seen to be a genuine property of the 
DMFT. The critical exponent 7 = 1 is in accordance 
with the mean-field nature of the DMFT, which neglects 
short-range spatial correlations between the electrons. 

In the similar manner one can show that at T = 0, 
where the transition becomes a quantum phase transi- 
tion, the static spin susceptibility diverges as the control 
parameter U approaches Uc(T ~ 0) from below 



X 



DMFT 



1 



(33) 



UciT = 0)-U' 
with the mean field exponent 7=1. 

While the dynamics, i.e., the effect of temporal cor- 
relations, is found to be very important in determining 
non-universal quantities such as the critical temperature 
or the Curie constant, we conclude that it apparently 
does not affect the universal scaling properties of the 
paramagnetic-to-ferromagnetic phase transition within 
the DMFT. 
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